Novel insights into genetic susceptibility for colorectal cancer from transcriptome-wide association and functional investigation.

BACKGROUND: Transcriptome-wide association studies have been successful in identifying candidate susceptibility genes for colorectal cancer (CRC). To strengthen susceptibility gene discovery, we conducted a large transcriptome-wide association study and an alternative splicing transcriptome-wide association study in CRC using improved genetic prediction models and performed in-depth functional investigations. METHODS: We analyzed RNA-sequencing data from normal colon tissues and genotype data from 423 European descendants to build genetic prediction models of gene expression and alternative splicing and evaluated model performance using independent RNA-sequencing data from normal colon tissues of the Genotype-Tissue Expression Project. We applied the verified models to genome-wide association studies (GWAS) summary statistics among 58 131 CRC cases and 67 347 controls of European ancestry to evaluate associations of genetically predicted gene expression and alternative splicing with CRC risk. We performed in vitro functional assays for 3 selected genes in multiple CRC cell lines. RESULTS: We identified 57 putative CRC susceptibility genes, which included the 48 genes from transcriptome-wide association studies and 15 genes from splicing transcriptome-wide association studies, at a Bonferroni-corrected P value less than .05. Of these, 16 genes were not previously implicated in CRC susceptibility, including a gene PDE7B (6q23.3) at locus previously not reported by CRC GWAS. Gene knockdown experiments confirmed the oncogenic roles for 2 unreported genes, TRPS1 and METRNL, and a recently reported gene, C14orf166. CONCLUSION: This study discovered new putative susceptibility genes of CRC and provided novel insights into the biological mechanisms underlying CRC development.

Genome-Wide Gene-Environment Interaction Analyses to Understand the Relationship between Red Meat and Processed Meat Intake and Colorectal Cancer Risk.

BACKGROUND: High red meat and/or processed meat consumption are established colorectal cancer risk factors. We conducted a genome-wide gene-environment (GxE) interaction analysis to identify genetic variants that may modify these associations. METHODS: A pooled sample of 29,842 colorectal cancer cases and 39,635 controls of European ancestry from 27 studies were included. Quantiles for red meat and processed meat intake were constructed from harmonized questionnaire data. Genotyping arrays were imputed to the Haplotype Reference Consortium. Two-step EDGE and joint tests of GxE interaction were utilized in our genome-wide scan. RESULTS: Meta-analyses confirmed positive associations between increased consumption of red meat and processed meat with colorectal cancer risk [per quartile red meat OR = 1.30; 95% confidence interval (CI) = 1.21-1.41; processed meat OR = 1.40; 95% CI = 1.20-1.63]. Two significant genome-wide GxE interactions for red meat consumption were found. Joint GxE tests revealed the rs4871179 SNP in chromosome 8 (downstream of HAS2); greater than median of consumption ORs = 1.38 (95% CI = 1.29-1.46), 1.20 (95% CI = 1.12-1.27), and 1.07 (95% CI = 0.95-1.19) for CC, CG, and GG, respectively. The two-step EDGE method identified the rs35352860 SNP in chromosome 18 (SMAD7 intron); greater than median of consumption ORs = 1.18 (95% CI = 1.11-1.24), 1.35 (95% CI = 1.26-1.44), and 1.46 (95% CI = 1.26-1.69) for CC, CT, and TT, respectively. CONCLUSIONS: We propose two novel biomarkers that support the role of meat consumption with an increased risk of colorectal cancer. IMPACT: The reported GxE interactions may explain the increased risk of colorectal cancer in certain population subgroups.

Identification of intergenerational epigenetic inheritance by whole genome DNA methylation analysis in trios.

Genome-wide association studies have identified thousands of loci associated with common diseases and traits. However, a large fraction of heritability remains unexplained. Epigenetic modifications, such as the observed in DNA methylation have been proposed as a mechanism of intergenerational inheritance. To investigate the potential contribution of DNA methylation to the missing heritability, we analysed the methylomes of four healthy trios (two parents and one offspring) using whole genome bisulphite sequencing. Of the 1.5 million CpGs (19%) with over 20% variability between parents in at least one family and compatible with a Mendelian inheritance pattern, only 3488 CpGs (0.2%) lacked correlation with any SNP in the genome, marking them as potential sites for intergenerational epigenetic inheritance. These markers were distributed genome-wide, with some preference to be located in promoters. They displayed a bimodal distribution, being either fully methylated or unmethylated, and were often found at the boundaries of genomic regions with high/low GC content. This analysis provides a starting point for future investigations into the missing heritability of simple and complex traits.

Corrigendum: CARB-ES-19 multicenter study of carbapenemase-producing Klebsiella pneumoniae and Escherichia coli from all Spanish provinces reveals interregional spread of high-risk clones such as ST307/OXA-48 and ST512/KPC-3.

[This corrects the article DOI: 10.3389/fmicb.2022.918362.].

Heterozygote advantage at HLA class I and II loci and reduced risk of colorectal cancer.

Objective: Reduced diversity at Human Leukocyte Antigen (HLA) loci may adversely affect the host's ability to recognize tumor neoantigens and subsequently increase disease burden. We hypothesized that increased heterozygosity at HLA loci is associated with a reduced risk of developing colorectal cancer (CRC). Methods: We imputed HLA class I and II four-digit alleles using genotype data from a population-based study of 5,406 cases and 4,635 controls from the Molecular Epidemiology of Colorectal Cancer Study (MECC). Heterozygosity at each HLA locus and the number of heterozygous genotypes at HLA class -I (A, B, and C) and HLA class -II loci (DQB1, DRB1, and DPB1) were quantified. Logistic regression analysis was used to estimate the risk of CRC associated with HLA heterozygosity. Individuals with homozygous genotypes for all loci served as the reference category, and the analyses were adjusted for sex, age, genotyping platform, and ancestry. Further, we investigated associations between HLA diversity and tumor-associated T cell repertoire features, as measured by tumor infiltrating lymphocytes (TILs; N=2,839) and immunosequencing (N=2,357). Results: Individuals with all heterozygous genotypes at all three class I genes had a reduced odds of CRC (OR: 0.74; 95% CI: 0.56-0.97, p= 0.031). A similar association was observed for class II loci, with an OR of 0.75 (95% CI: 0.60-0.95, p= 0.016). For class-I and class-II combined, individuals with all heterozygous genotypes had significantly lower odds of developing CRC (OR: 0.66, 95% CI: 0.49-0.87, p= 0.004) than those with 0 or one heterozygous genotype. HLA class I and/or II diversity was associated with higher T cell receptor (TCR) abundance and lower TCR clonality, but results were not statistically significant. Conclusion: Our findings support a heterozygote advantage for the HLA class-I and -II loci, indicating an important role for HLA genetic variability in the etiology of CRC.

Improving Species Level-taxonomic Assignment from 16S rRNA Sequencing Technologies.

Analysis of the bacterial community from a 16S rRNA gene sequencing technologies requires comparing the reads to a reference database. The challenging task involved in annotation relies on the currently available tools and 16S rRNA databases: SILVA, Greengenes and RDP. A successful annotation depends on the quality of the database. For instance, Greengenes and RDP have not been updated since 2013 and 2016, respectively. In addition, the nature of 16S sequencing technologies (short reads) focuses mainly on the V3-V4 hypervariable region sequencing and hinders the species assignment, in contrast to whole shotgun metagenome sequencing. Here, we combine the results of three standard protocols for 16S rRNA amplicon annotation that utilize homology-based methods, and we propose a new re-annotation strategy to enlarge the percentage of amplicon sequence variants (ASV) classified up to the species level. Following the pattern (reference) method: DADA2 pipeline and SILVA v.138.1 reference database classification (Basic Protocol 1), our method maps the ASV sequences to custom nucleotide BLAST with the SILVA v.138.1 (Basic Protocol 2), and to the 16S database of Bacteria and Archaea of NCBI RefSeq Targeted Loci Project databases (Basic Protocol 3). This new re-annotation workflow was tested in 16S rRNA amplicon data from 156 human fecal samples. The proposed new strategy achieved an increase of nearly eight times the proportion of ASV classified at the species level in contrast to the reference method for the database used in the present research. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Sample inference and taxonomic profiling through DADA2 algorithm. Basic Protocol 2: Custom BLASTN database creation and ASV taxonomical assignment. Basic Protocol 3: ASV taxonomical assignment using NCBI RefSeq Targeted Loci Project database. Basic Protocol 4: Definitive selection of lineages among the three methods.

Identification of a Twelve-microRNA Signature with Prognostic Value in Stage II Microsatellite Stable Colon Cancer.

We aimed to identify and validate a set of miRNAs that could serve as a prognostic signature useful to determine the recurrence risk for patients with COAD. Small RNAs from tumors of 100 stage II, untreated, MSS colon cancer patients were sequenced for the discovery step. For this purpose, we built an miRNA score using an elastic net Cox regression model based on the disease-free survival status. Patients were grouped into high or low recurrence risk categories based on the median value of the score. We then validated these results in an independent sample of stage II microsatellite stable tumor tissues, with a hazard ratio of 3.24, (CI95% = 1.05-10.0) and a 10-year area under the receiver operating characteristic curve of 0.67. Functional analysis of the miRNAs present in the signature identified key pathways in cancer progression. In conclusion, the proposed signature of 12 miRNAs can contribute to improving the prediction of disease relapse in patients with stage II MSS colorectal cancer, and might be useful in deciding which patients may benefit from adjuvant chemotherapy.

TALKIEN: crossTALK IntEraction Network. A web-based tool for deciphering molecular communication through ligand-receptor interactions.

Molecular crosstalk, the dialogue between different cell types, is attracting more attention in cancer research. On the one hand, the communication between tumor and non-tumor cells in the microenvironment or between different tumor clones has influential consequences for the progression and spread of tumors and response to treatment. On the other hand, novel techniques such as single-cell sequencing or spatial transcriptomics provide detailed information that needs to be interpreted. TALKIEN: crossTALK IntEraction Network is a simple and intuitive online R/shiny application to visualize molecular crosstalk information through the construction and analysis of a protein-protein interaction network. Taking two or more lists of genes or proteins as input, which are representative of cell lineages, TALKIEN extracts information about ligand-receptor interactions, builds a network and analyzes it using systems biology techniques such as centrality measures and component analysis, among others. Moreover, it expands the network displaying pathways downstream receptors. The application allows users to select different graphical layouts, performs functional analysis and gives information about drugs targeting receptors. In conclusion, TALKIEN allows users to detect ligand-receptor interactions generating new in silico predictions of cell-cell communication thus providing a translational rationale for future experiments. It is freely available at https://www.odap-ico.org/talkien.

Recommendations of the Spanish Antibiogram Committee (COESANT) for in vitro susceptibility testing of antimicrobial agents by disk diffusion.

Disk diffusion is a well standardized method that provides reliable categorical results to guide antimicrobial therapy in numerous types of infections. Based on the guidelines of the European Committee on Antimicrobial Susceptibility Testing (EUCAST), which are widely implemented in Spain, the Spanish Antibiogram Committee (COESANT) has drawn up recommendations for antimicrobial selection by the disk diffusion technique, including selective reporting and its use for the detection of resistance mechanisms. Factors affecting disk diffusion results, along with advantages and shortcomings of the method, are also discussed.

Spanish Antibiogram Committee (COESANT) recommendations for cumulative antibiogram reports.

The Spanish Antibiogram Committee (Comité Español del Antibiograma, COESANT) presents in this document a series of recommendations intending to unify how cumulative antibiogram reports must be made in Clinical Microbiology Spanish laboratories. This article is based on the information included in the Clinical Microbiology Procedure No. 51, «Preparation of cumulative reports on antimicrobial susceptibility¼ of the Spanish Society of Infectious Diseases and Clinical Microbiology (SEIMC), published in 2014. The recommendations also include the modifications in the definition of clinical interpretive categories recently published by the European Committee on Antimicrobial Susceptibility Testing (EUCAST) in 2019. Its final objective is to establish a homogeneous way of preparing these summaries to compare results from different centers or aggregate the information from these in order to carry out an adequate local or even national surveillance regarding the evolution of antimicrobial susceptibility.

Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries.

Colorectal cancer (CRC) is a leading cause of mortality worldwide. We conducted a genome-wide association study meta-analysis of 100,204 CRC cases and 154,587 controls of European and east Asian ancestry, identifying 205 independent risk associations, of which 50 were unreported. We performed integrative genomic, transcriptomic and methylomic analyses across large bowel mucosa and other tissues. Transcriptome- and methylome-wide association studies revealed an additional 53 risk associations. We identified 155 high-confidence effector genes functionally linked to CRC risk, many of which had no previously established role in CRC. These have multiple different functions and specifically indicate that variation in normal colorectal homeostasis, proliferation, cell adhesion, migration, immunity and microbial interactions determines CRC risk. Crosstissue analyses indicated that over a third of effector genes most probably act outside the colonic mucosa. Our findings provide insights into colorectal oncogenesis and highlight potential targets across tissues for new CRC treatment and chemoprevention strategies.

COLONOMICS - integrative omics data of one hundred paired normal-tumoral samples from colon cancer patients.

Colonomics is a multi-omics dataset that includes 250 samples: 50 samples from healthy colon mucosa donors and 100 paired samples from colon cancer patients (tumor/adjacent). From these samples, Colonomics project includes data from genotyping, DNA methylation, gene expression, whole exome sequencing and micro-RNAs (miRNAs) expression. It also includes data from copy number variation (CNV) from tumoral samples. In addition, clinical data from all these samples is available. The aims of the project were to explore and integrate these datasets to describe colon cancer at molecular level and to compare normal and tumoral tissues. Also, to improve screening by finding biomarkers for the diagnosis and prognosis of colon cancer. This project has its own website including four browsers allowing users to explore Colonomics datasets. Since generated data could be reuse for the scientific community for exploratory or validation purposes, here we describe omics datasets included in the Colonomics project as well as results from multi-omics layers integration.

Multimorbidity and chronic co-prescription networks and potential interactions in adult patients with epilepsy: MorbiNet study.

OBJECTIVE: We constructed epilepsy multimorbidity networks to study associations with chronic conditions, and co-prescriptions and drug-disease networks to assess potential interactions. We conducted a population-based study in Catalonia, Spain, with electronic files of 3,135,948 adult patients with multimorbidity, 32,625 of them with epilepsy (active diagnosis any time during 2006-2017). We constructed epilepsy comorbidity networks using logistic regression models from odds ratio estimates adjusted by age, sex, and comorbidities with R software and generated trajectories to study the progression of epilepsy. We constructed drug-disease and co-prescription networks using mixed models with repeated measures adjusting by age, sex, and period with chronic prescription invoiced data. Comorbidity more frequently preceding epilepsy included cerebrovascular accident (OR: 3.59), congenital anomalies (2.18), and multiple sclerosis (1.33); and following epilepsy: dementia (1.91), personality disorder (1.59), alcohol abuse (1.22), and Parkinson (1.21). Mental retardation (13.08), neurological cancer (8.49), benign neoplasm (4.69), infections (3.14), and psychosis (1.58) might precede or not epilepsy. A common progression was to schizophrenia, dementia, and other neurological diseases (mainly cerebral palsy and other degenerative diseases of nervous system). Co-prescription associations with major-moderate potential interactions were 54% for carbamazepine, 61% phenytoin, 53% phenobarbital, and 32% valproate. Major potential interactions were with antipsychotic, anxiolytic, opioid, cardiovascular, and other anti-seizure medications (ASMs). The most frequent comorbidities of epilepsy were congenital, cerebrovascular, and neurological and psychiatric conditions. High comorbidity and co-prescription with potential interactions can increase the complexity of care of patients with epilepsy.

Meta-Analysis and Validation of a Colorectal Cancer Risk Prediction Model Using Deep Sequenced Fecal Metagenomes.

The gut microbiome is a potential modifiable risk factor for colorectal cancer (CRC). We re-analyzed all eight previously published stool sequencing data and conducted an MWAS meta-analysis. We used cross-validated LASSO predictive models to identify a microbiome signature for predicting the risk of CRC and precancerous lesions. These models were validated in a new study, Colorectal Cancer Screening (COLSCREEN), including 156 participants that were recruited in a CRC screening context. The MWAS meta-analysis identified 95 bacterial species that were statistically significantly associated with CRC (FDR < 0.05). The LASSO CRC predictive model obtained an area under the receiver operating characteristic curve (aROC) of 0.81 (95%CI: 0.78−0.83) and the validation in the COLSCREEN dataset was 0.75 (95%CI: 0.66−0.84). This model selected a total of 32 species. The aROC of this CRC-trained model to predict precancerous lesions was 0.52 (95%CI: 0.41−0.63). We have identified a signature of 32 bacterial species that have a good predictive accuracy to identify CRC but not precancerous lesions, suggesting that the identified microbes that were enriched or depleted in CRC are merely a consequence of the tumor. Further studies should focus on CRC as well as precancerous lesions with the intent to implement a microbiome signature in CRC screening programs.

Host-associated variability of the cdtABC operon, coding for the cytolethal distending toxin, in Campylobacter jejuni.

Campylobacter, a major cause of food-borne gastroenteritis worldwide, colonize the gastrointestinal tract of a wide range of animals, being birds the main reservoir. The mechanisms involved in the interaction of Campylobacter with the different hosts are poorly understood. The cytolethal distending toxin, encoded in the cdtABC operon, is considered a pivotal virulence factor during human infection. Differences in the prevalence of cdtABC genes in Campylobacter isolates from three distinct origins (wild birds, broiler chickens and humans) prompted us to further characterize their allelic variability. The sequence of cdtABC is highly conserved among broiler and human isolates. A high diversity of cdtABC alleles was found among wild bird isolates, including several alleles that do not produce any functional CDT. These results suggest that specific variants of the cdtABC operon might define the host range of specific Campylobacter jejuni isolates. Moreover, our data indicate that PCR methodology is inaccurate to characterize the prevalence of the cdt genes, since negative PCR detection can be the result of divergences in the sequence used for primer design rather than indicating the absence of a specific gene.

CARB-ES-19 Multicenter Study of Carbapenemase-Producing Klebsiella pneumoniae and Escherichia coli From All Spanish Provinces Reveals Interregional Spread of High-Risk Clones Such as ST307/OXA-48 and ST512/KPC-3.

Objectives: CARB-ES-19 is a comprehensive, multicenter, nationwide study integrating whole-genome sequencing (WGS) in the surveillance of carbapenemase-producing K. pneumoniae (CP-Kpn) and E. coli (CP-Eco) to determine their incidence, geographical distribution, phylogeny, and resistance mechanisms in Spain. Methods: In total, 71 hospitals, representing all 50 Spanish provinces, collected the first 10 isolates per hospital (February to May 2019); CPE isolates were first identified according to EUCAST (meropenem MIC > 0.12 mg/L with immunochromatography, colorimetric tests, carbapenem inactivation, or carbapenem hydrolysis with MALDI-TOF). Prevalence and incidence were calculated according to population denominators. Antibiotic susceptibility testing was performed using the microdilution method (EUCAST). All 403 isolates collected were sequenced for high-resolution single-nucleotide polymorphism (SNP) typing, core genome multilocus sequence typing (cgMLST), and resistome analysis. Results: In total, 377 (93.5%) CP-Kpn and 26 (6.5%) CP-Eco isolates were collected from 62 (87.3%) hospitals in 46 (92%) provinces. CP-Kpn was more prevalent in the blood (5.8%, 50/853) than in the urine (1.4%, 201/14,464). The cumulative incidence for both CP-Kpn and CP-Eco was 0.05 per 100 admitted patients. The main carbapenemase genes identified in CP-Kpn were bla OXA-48 (263/377), bla KPC-3 (62/377), bla VIM-1 (28/377), and bla NDM-1 (12/377). All isolates were susceptible to at least two antibiotics. Interregional dissemination of eight high-risk CP-Kpn clones was detected, mainly ST307/OXA-48 (16.4%), ST11/OXA-48 (16.4%), and ST512-ST258/KPC (13.8%). ST512/KPC and ST15/OXA-48 were the most frequent bacteremia-causative clones. The average number of acquired resistance genes was higher in CP-Kpn (7.9) than in CP-Eco (5.5). Conclusion: This study serves as a first step toward WGS integration in the surveillance of carbapenemase-producing Enterobacterales in Spain. We detected important epidemiological changes, including increased CP-Kpn and CP-Eco prevalence and incidence compared to previous studies, wide interregional dissemination, and increased dissemination of high-risk clones, such as ST307/OXA-48 and ST512/KPC-3.